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Abstract 

We express the two-massless-flavor Gell-Mann-Oakes-Renner ratio in terms of 
low-energy tttt observables including the 0(p 6 ) double chiral logarithms of general- 
ized chiral perturbation theory. Their contribution is sizeable and tends to com- 
pensate the one from the single chiral logarithms. However it is not large enough 
to spoil the convergence of the chiral expansion. As a signal of reduced theoretical 
uncertainty, we find that the scale dependence from the one-loop single logarithms 
is almost completely canceled by the one from the two-loop double logarithms. 
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1. The dynamics of chiral symmetry breakdown (%SB) in the Standard Model remains 
an open issue both from the experimental and theoretical standpoint. At very low energy the 
problem concerns the understanding of the mechanism of for QCD. In this framework 
it is generally assumed that chiral symmetry is broken through the formation of a quark- 
antiquark condensate. A theoretically consistent alternative is that the symmetry is broken 
by higher dimensional order parameters, analogously to what happens in antiferromagnetic 
spin systems for the rotational symmetry ||TJ. Deeper insight into this problem can be gained 
by considering the theory in a Euclidean box, with (anti)periodic boundary conditions. The 
two alternatives are then distinguished by a qualitatively different behavior of the smallest 
eigenvalues of the Dirac operator in the thermodynamical limit. In fact it has recently been 
realized that a prominent role is played by the number of massless flavors Nf of the theory 
P, H|]. In general, as Nf increases, the order parameters dominated by the infrared end of 
the spectrum of the Euclidean Dirac operator experience a paramagnetic suppression which, 
for Nf large enough, will eventually restore chiral symmetry. Since only the SU(2) and SU(3) 
chiral limits are reasonable approximations to the real world, the question to be addressed 
is to which extent Nf = 2, 3 is close to a critical point in which (qq) disappears and chiral 
symmetry is eventually broken by higher-dimensional order parameters. As was pointed out in 
Refs. [f|, |J, it is possible to submit the standard hypothesis of a large condensate to experimental 
verification, through the measurement of low-energy phase-shifts. The S-wave tttt interaction 
at low energy is very sensitive to the two-massless-flavor chiral condensate: it gets stronger as 
the condensate decreases. Standard chiral perturbation theory (ChPT) || pushed to two- 
loop accuracy, together with the use of constraints from Roy equations |7J, provides a very 
sharp prediction for the S-wave scattering lengths ||. The situation has evolved recently, 
since new preliminary experimental data provided by the E865 experiment on i^ e4 decays have 
become available ||. In fact a first analysis of the above mentioned data in the standard 
context indicates that deviations from this scenario are marginal [fLOfl . On the other hand, 
generalized chiral perturbation theory (GChPT), which relaxes from the beginning the standard 
assumption, cannot produce any prediction for the low-energy Tin observables, but only relates 
these observables to the magnitude of the quark condensate. It is the purpose of the present 
work to establish the relationship between low-energy tctt observables and the two-massless- 
flavors quark condensate in the language of GChPT beyond the 0(p A ) level, i.e. including 
the 0(p 6 ) double chiral logarithms, which are among the potentially most dangerous two-loop 
contributions. 

2. Relaxing the standard assumption of a dominant condensate implies that the Gell- 
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Mann-Oakes-Renner ratioQ, 



is actually a free parameter. If its value were not close to 1, then the quadratic term in 
quark masses would be as important as the condensate term in the chiral expansion of the 
squared pion mass. Thus the counting rule has to be modified accordingly: both the light 
quark mass in (m u = rrid = rh) and the quark condensate are quantities of chiral order 0(p). 
The effective Lagrangian is organized as an infinite sum C = + + + . . . where the 
0(p d ) term contains terms with k derivatives, I powers of quark masses and n powers of 
the quark condensate, such that k + l + n = d. It contains, at each finite chiral order, additional 
operators with respect to the corresponding standard Lagrangian. For instance, the leading 
0{p 2 ) Lagrangian reads, 

£ (2) = ^-{(D^D»U) + 2B(U^x + X 1 U)+A((tfx) 2 + (x f u) 2 ) 

+Z P (U^ X - X ] U) 2 + h (x^x) + K (det X + det X ] ) } , (2) 

where U G SU(2) collects the Goldstone bosons' degrees of freedom and x denotes the scalar- 
pseudoscalar external source, to be expanded around the quark mass matrix x = ^1 + • • •• 
Notice that, compared to standard notations, a factor 2B has been removed from the definition 
of Xi f° r consistency with the new chiral counting: in fact B is related to the quark condensate, 
F 2 B = — lim rfl _> (gg), and it has to be considered, formally, as a small parameter, B ~ x ~ 
0(p). One- loop diagram with insertions from the Lagrangian (fj) are of chiral order 0(p 4 ). 
They are divergent and require the renormalization of the low-energy constants (I.e. a). In the 
minimal subtraction scheme of dimensional regularization, 

fi d ~ A r const 

const = const 7 * + — — , cons | , (3) 

where \x is the renormalization scale. The /3-function coefficients r const for the leading l.e.c.'s 
are 

F 2 T A = 3B 2 , F 2 T ZP = ~B 2 , F 2 T h , = QB 2 , (4) 

whereas F 2 , B and h do not get renormalized. The renormalized constants are renormalization 
scale dependent, according to 

;U-^-const r = -— ^jr const . (5) 



1 In Eq. (|l|), F%M% should be considered simply as physical units for measuring the QCD renormalization 
group invariant combination (m u + md){qq). 
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Due to the chiral counting of the condensate parameter B, the divergences arising from the 
renormalization of the l.e.c.'s are of chiral order p 4 , as they should, in order to absorb the ones 
arising from the loops. The renormalization of the leading l.e.c.'s (|j) is not sufficient to absorb 
all the one-loop 0(p i ) divergences generated by C^ 2 \ which is a consequence of the fact that 
the Lagrangian (0) is non-renormalizable. One has to include also the subleading 0(p 3 ) and 
0(p 4 ) Lagrangians and renormalize the relative l.e.c.'s, according to Eq. (Q). The /3-function 
coefficients T C (3) for the l.e.c.'s of will be proportional to B, in order that the divergences 
be of order p 4 , while those for the l.e.c.'s of will be independent of B, 

r £(2) ~ o(b 2 ) ~ o( P 2 ), r £(3) ~ o(B) ~ o( P ), r cW ~ 0(1). (6) 

We refer to existing literature for the complete list of operators together with their renormal- 
ization up to 0(p 4 ) p"T|j. 



3. The proliferation of chirally invariant operators has no bearing on the form of the 
mi scattering amplitude up to two-loop level: the latter only depends on six "observable" 
parameters, independently of the chiral counting. This is a consequence of unitarity, analyticity 
and crossing symmetry, and of the chiral suppression of partial waves greater than two 0. 
Neglecting 0(p 8 ), i.e. at two- loop level, the invariant amplitude for txtx scattering has been 
given, using dispersive methods, in Ref. [0 in terms of the parameters a, (3, Ai, . . . , A4, 

A(s\t, u) = A KMSF (s\t, u; a, (3; X h A 2 , A 3 , A 4 ) + 0(p 8 ). (7) 

The parameter a represents, at leading order, the amplitude at the symmetric point s = 
t = u = 4/3M 2 , and j3 its slope. These two parameters can be viewed as the subtraction 
constants in a Roy-like dispersive representation of the amplitude which fixes, confronted to 
existing experimental data in the medium and high energy region, the remaining 4 parameters 
Ai, . . . , A4 A brute force calculation in ChPT allows one to find the "chiral anatomy" of 



these parameters, i.e. their expressions in terms of quark masses and l.e.c.'s. In particular 
we are interested in finding their dependence on the l.e.c. B, related to the quark condensate. 
At tree level only the parameter a depends on the condensate: its value varies from 1 to 4 
if (qq) decreases from its standard value down to zero, while (3 is equal to 1 and the Aj's 
vanish. Going to one-loop precision already involves a conspicuous set of l.e.c. characteristics 
of GChPT, besides the contribution of the chiral logarithms. Explicit expressions for the 1-loop 



parameters a, (3, Ai and A2 can be found in Refs. |lT|. After eliminating as many unknown 
constants as possible in favor of physical quantities, one can invert these expressions for the 
quark condensate, and write the GOR ratio in terms of the combination a + 2/3, modulo 



remaining unknown constants, renormalized at a scale /1 [|14 

2mBF 2 a + 2(3 F 2 

FlMl 3 + F1MI 
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+ 



Aa r 2 + 8 - 1 j a 3 + 8b r 2 + lQc[ 



rh 2 



64 

7T 

+ 288^ (a + 2/?) t 24 - 11 (« + 2 / 3 )] 



+ 



6 + jj(a + 2/3)-H( a + 2 /3) 2 



M 2 . M; 
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32tt 2 F 2 ° /i 2 



log • (8) 



The role of a and /3 in this equation is that of observable quantities to be extracted from 
the data by fitting the two- loop formula (|7|), with the A, fixed by the Roy equations. The 
standard prediction for a and (3 is very close to 1 |T5| . On the other hand the contribution 
from the l.e.c.'s figuring in Eq. @ is unknown, at present. The pj's come from the 0(p 3 ) 
£(0,3) Lagrangian, with no derivatives and three powers of the scalar-pseudoscalar source, the 
constants Oj, 6, and Cj come from the 0(p 4 ) £(2,2) Lagrangian and the constants e;, /j from 
the 0(p 4 ) £(0,4) Lagrangian, with analogous notations. In principle these l.e.c.'s should be 
determined from independent phenomenological information but, due to the large number of 
independent operators, this program seems hopeless. However we do know the scale dependence 
of all these l.e.c.'s. One possibility is thus to treat them as randomly distributed around zero 
with an error assigned to them according to naive dimensional analysis estimates, 

Pi ~ ±y-, a i ,b i ,c i ,e i ,f i ~±-^- J (9) 
Ah A h 

where the hadronic scale Ah will be fixed to 1 GeV. The light quark mass rh is written as 
rh = m s /r, where the quark mass ratio r is strongly correlated to the two-flavor GOR ratio (see 
e.g. Fig. 1 of Ref. ||) and m s is taken conservatively to be ~ 200 MeV. As a consequence of 
the smallness of the light quark mass, these contributions will in practice be very small. On the 
contrary, the logarithmic terms in Eq. (|8|) will be quite important, specially in the region where 
a + 2/3 is largeQ. A further source of uncertainty is therefore the scale \x inside the logarithm, 
i.e. the scale at which the estimates (H) are supposed to hold. We will take \i = M^±250 MeV. 
The resulting curve is shown as the dashed line in Fig. [I], the band representing the theoretical 
uncertainty obtained as explained above. The main contribution to the error band is represented 
by the scale variation inside the logarithm. Due to the large contribution of the latter it 
becomes necessary to test the convergence of the chiral series, by studying the importance of 
higher orders. In the following we will compute the 0(p 6 ) double logarithmic corrections to 
this result. 



2 Notice that the logarithm is absent in the standard case, in which a = j3 = 1 in higher orders. Indeed a 
and (3 are scale-independent quantities. Since all unknown l.e.c.'s in Eq. are at least 0(p 6 ) in the standard 
counting, there is no possibility for a standard chiral logarithm at order p 4 . 
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4. Double chiral logarithms are among the potentially most dangerous contributions at 
order 0(p e ) ||16|| , because of the smallness of the pion mass. As first pointed out in Ref. [17 



(see also Ref. |T8|) they can be obtained from a 1-loop calculation, using the fact that, in the 
renormalization procedure, non-local divergences must cancel. We set the space-time dimension 
d = 4 + uj to regulate the theory. The coupling constant F 2 , which multiplies the kinetic term 
for the pions, has dimension [F 2 ] = d — 2. This constant appears as a common factor in the 
definition of the generalized chiral Lagrangian. Therefore, a generic I.e. a, with k derivatives 
and I powers of the scalar-pseudoscalar source, c^k,i), has mass dimension 2 — k — I, and does 
not depend on to. We thus replace F 2 with fi^F 2 , making appear explicitly the scale parameter 
H brought in by the regularization procedure. Since each loop involves a factor F~ 2 , the chiral 
expansion of a generic amplitude A up to two loops, apart from an overall dimensional factor, 
takes the form, 

^Ve^ (1) (^)^ 1oop + ( —Ty:p?\c a )fr oop , m 



(1 2) 

where P-; ' are polynomials in the quark mass and l.e.c.'s, generically denoted by c a , and gi and 
fi are loop-functions of the kinematical variables, expressed in terms of dimensionless quantities 
with the appropriate inverse powers of M n . The loop functions are singular as ui — > 0, 

fi = k2 + ki +f ■ +o{u>), ft = ^ + Sil o + o(4 (ii) 

UJ A UJ UJ 

While the one-loop divergences are canceled by the renormalization of the l.e.c.'s of + + 
£ (4) , appearing in Aree, 

c a = - + c r a + 0(uj), (12) 

UJ 

the two-loop divergences require the introduction of higher order local counterterms. We have 
defined, compared to the notations of Eq. (|3|), 

167T 2 

which is actually independent of uj, due to the fact that the /3-function coefficients T a are all 
proportional to F~ 2 (cfr. Eq. (f|) and Ref. [|ll|] for the and l.e.c.'s). Notice also that 
the divergent term in Eq. ( |T2"D always arises at chiral order p A in the generalized counting, as 
already mentioned, even if the constant itself is of lower order. Therefore, since the 1-loop 
polynomials are already of order p 4 , we never have to deal, up to and including order p e , 
with products of the type 5i5jg 1 ~ loop : the polynomials P^ after the substitution (|T2"D, can be 
rewritten as 

pPM = H 1} K) + ^K)^. (14) 



< I rZ 
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For the same reason, since the 2-loop functions are of order p 6 , the l.e.c.'s inside the polynomials 

(2) 

can be simply replaced by the corresponding renormalized ones, 



P 



(2), 



P, 



( 2 ) ( 



(15) 



With the replacements ( |TTD and (|12|) the u>-dependence of Eq. ( |10|) is explicit. General theorems 
of renormalization theory require the residues of the poles in uo to be polynomials in external 
momenta and masses. Thus the non-local divergences, of the type (\ogM n )/u must cancel in 
the final result. This condition amounts to a relation between the coefficients and g^x, 
namely 



< dc b V 



(16) 



This is the reason why the double logarithms can be obtained by a 1-loop calculation: they 
only enter through the residues g^i of the 1-loop functions. After renormalization the generic 
amplitude A becomes, 



A = A r ree + E{ai 



(MS 




+9i,o 



+E{^ (2) «; 



+ E^(^log(^ 



Ao + 2/^logf^ 



;i7) 



Thus we have to compute, besides the 1-loop graphs with the 0{p 2 ) Lagrangian, all the 1- 
loop graphs with 1 insertion of operators from the 0(p 4 ) Lagrangian and up to 2 insertions of 
operators from the 0(p 3 ) Lagrangian. At the end, each occurrence of log(M 7r //i) must be 
replaced by the combination 

1 (M T 
7*- log -f 




1 




and all contributions beyond 0(p 4 ), but the pure double logarithms, must be discarded. Obvi- 
ously this is only part of the 0(p 6 ) contribution. Let us examine in more detail what we are 
considering and what we are neglecting by this procedure. The double logarithms are of chiral 
order p 6 and always arise paired with single logarithms multiplied by a I.e. a, as in Eq. 



They are always of equal or higher order, compared to the accompanying single logarithms: the 
latters, in the generalized counting, can be of order p 4 , p 5 or p 6 . We are keeping obviously the 
p 4 single logarithms, but neglecting the p 5 and p 6 . In any case, we could only consider them 
as error sources, since most of the and l.e.c.'s are unknown experimentally. The fact 
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that, at 0(p 4 ), the error due to the unknown constants is much smaller than the one coming 
from the variation of the scale in the logarithms, is one argument in favor of keeping only the 
pure double logarithms. On the other hand, at least in the standard case, the pure double 
logarithmic contribution to a is by a factor 10 larger than the contributions of the type l^L 



Uj. We are only interested in the relationship between the GOR ratio and the parameters a 



and f3. Since the parameter a is the most correlated to the quark condensate, we can expect 
that the latter also (and therefore the GOR ratio) be dominated, at 0(p 6 ), by the pure double 
logarithms. 



5. The first step for computing the amplitude is the calculation of the axial- axial 
two-point function, from which we can extract F n and M n . We display the result for these two 
quantities, where all l.e.c.'s, here and in the following, are renormalized at a scale p, but we 
drop the superscript r for simplicity: 



F 2 

—M 2 
pi * 



2Bm + AArh 2 + (9p 1 + p 2 + 20p 4 + 2p 5 ) m 3 



+ (16 ei + 4e 2 + 32/! + 40/ 2 + 8/ 3 + 96/ 4 ) m 
+4a 3 M 2 m 2 - (3M, 



32tt 2 F 2 



2 ' 20 Am 2 ) log ^$ 



+ 



33 65 Am 2 ( Am 2 \ 

1 + 2£ {2) m + (2ai + a 2 + 4a : 
Ml Ml 1MU 1 

8n 2 F 2 %2 ^ 2 FM 167T 2 



nr2 f M 2 , M 2S 
M log — 



26i - 26 2 ) m 2 




(19) 



(20) 



While the double chiral logarithm for F n is the same as in the standard case, for there are 
additional double logarithmic contributions which would be relegated by the standard counting 
at orders p 8 and p 10 . Whether these additional corrections are important or not, depends on 
the ratio Am 2 /Ml, i.e. on the deviation of the GOR ratio from 1. The tctt amplitude can be 
brought to the form of Eq. (|7|), explicitly displayed in Ref. []12~ |, with the following expressions 
for the parameters a, (3, Ai, . . . , A4 in terms of the l.e.c.'s 3 : 

F 2 



a 



F 2 M 2 



l.2Bm + IQAm 2 - AM 2 ^ 2) m + (81p x + p 2 + 164p 4 + 2p 



5 ) m 



-8M 2 (26i - 26 2 - a 3 - 4 Cl ) m 2 

+16 {6Aa 3 + 16ea + e 2 + 32^ + 34/ 2 + 2f 3 + 72/ 4 ) m 4 

A 2 m v 



Ml 
32n 2 F 2 



AM 2 + 20AAm 2 + 528 




i These results were formerly reported in Ref. pi 
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32vr 2 F 2 
533 



+M. 



+5808 



M 4 + 88Am 2 M* + 528A 2 m A 

Am 2 



18817 Am 2 61076 
72 + 30 Ml + 15 



M 2 



im 2 
~Mj 



Ml 



Ml 



(21) 



1 + 2£ (2) m - 4^ (2)2 m 2 + 2 (3a 2 + 2a 3 + 46i + 26 2 + 4ci) m 2 



4M 2 



32vr 2 F 2 
/151 



1 + 10 



400 




f + 1 



+ ( -^M 4 + — M 2 Am 2 + 420A 2 m 4 
V 36 3 



1 



Ai 



A, 




16tt 2 F 2 
130 Afn 2 ^ 



9 M 2 
80 im 2 ^ 



Ml 



M 2 

"I 7T 

16tt 2 F 2 ° g ~^~_ 



3 + 9 Ml 



Ml 



1Qti 2 FI 



log 



Ml 
V 2 



167T 



log 



Jig 
/z 2 



(22) 
(23) 
(24) 
(25) 



It is easy to check that these formulae, when restricted to the standard case, agree with the 
ones displayed in Ref.[I3| based on the complete two-loop calculation | J20| . 

Eliminating the constant A in favor of a and the £(2,1) coupling ^ in favor of (3, one can 
express the GOR ratio of Eq. (fj) as a function of the combination a + 2(3, 



x 



GOR 



X 



GOR 



1— loop 



+ 



121 
108 



a + 2pf 



4169, 5639, 

(a + 2/?)+ 1620 (a + 2/5) 

+ (M 2 m log M 2 ,M 2 m 




(26) 



where x^° R |i-ioop is the 0(p 4 ) result corresponding to Eq. (||). We are neglecting single- 
logarithmic and constant contributions which start at order p 5 in the generalized counting (cfr. 
discussion at the end of the previous section). The function (p6|) is shown in Fig. |] as the 
full line (lower curve), together with the theoretical uncertainty, represented by the shaded 
band and estimated varying the ChPT renormalization scale fi by ±250 MeV around the value 
fi = M v , with the unknown l.e.c.'s considered as error sources, as in Eq. ([§). In the lower band 
also, the error comes predominantly from the variation of the scale. 



6. As it is clear from Fig. [I], the contribution of the double logarithms to 



£2 is quite 

large, specially in the "extreme" generalized case of large a + 2/3, and tends to compensate 
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Figure 1: The GOR ratio at 1 loop including (continuous) and not including (dashed) the double 
logarithms. The scale is set to fi = M v ± 250 MeV. 

the 0(p A ) single logarithms. Still it is small enough in order to maintain the validity of the 
chiral expansion, the 0(p 6 ) weighting about a half of the 0(p A ) corrections. Moreover, the 
scale dependence is almost completely canceled by the double logarithms (cfr. the width of 
the shaded bands around the two curves). We stress that, at 0{p A ) the scale-dependence 
arises because of our ignorance about the l.e.c.'s. It reflects the uncertainty about the scale 
at which the estimates @ are supposed to hold. Indeed the effective theory, up to a given 
order of accuracy, defines a perfectly consistent theory, in the sense that the scale dependence 
from the chiral logarithms is compensated by the scale dependence of the l.e.c.'s [see Eq. (|5|)]. 
The non-renormalizability of the theory means that, as we want to increase the accuracy, we 
are forced to consider more and more l.e.c.'s, which implies a loss of predictive power of the 
theory. Therefore, if we knew the values of the l.e.c.'s at some scale, the complete 0(p A ) result 
would be scale- independent. On the contrary, an 0(p 6 ) calculation in the double logarithmic 
approximation introduces a scale dependence, which would be absorbed by the remaining 0(p 6 ) 
corrections, including further unknown l.e.c.'s. It is remarkable that these two distinct sources 
of scale-dependence almost completely cancel with each other. This is certainly a signal of 
reduced theoretical uncertainty, although one must admit that the cancellation is fortuitous. A 
complete two-loop calculation with the generalized Lagrangian, would not be very useful, due to 
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the flourishing of unknown l.e.c.'s. Double chiral logarithms are only part of this full two-loop 
calculation, but their contribution is known and unambiguous, apart from the scale-dependence. 
Moreover, there are reasons to believe that for some observables, namely the S-wave, isoscalar 
scattering length, or the parameter a, they constitute -at least in the standard case, where a 
complete 0(p 6 ) calculation is available- the bulk of the 0(p e ) corrections |16|]. We remind that 



the status of a+2/3 in this relationship is that of an observable quantity, extracted from the data 
using the full two-loop six-parametrical formula ([7]), with the constants \i(a, (3) determined e.g. 



from the recent solution of the Roy equations [51]. This procedure yields, through Eq. 



an experimental value for the two-flavor GOR ratio, with a theoretical uncertainty that we 
estimate to be twice the width of the lower band of Fig. [l]. It is worth recalling that a GChPT 
fit to the old Rosselet data |55[ gives a = 2.16 ± 0.86, (3 = 1.074 ± 0.053 |[5J, where the error 
was dominated by the statistics, and no correlations between the two constants where taken 
into account. The experimental uncertainty will certainly be reduced by the new data from the 
K e 4 experiment, E865 at Brookhaven M, with tenfold statistics compared to the Cern-Munich 
|22]| . A further increase in statistics is also expected in the planned NA48-2 K ei experiment 



one 



at CERN, which will provide a very welcome complementary set of data in terms of different 
systematic uncertainties [^3] . 
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